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ABSTRACT 

Over the last decade, evidence has mounted that the solar system's observed state can be favorably 
reproduced in the context of an instability-driven dynamical evolution model, such as the "Nice" 
model. To date, all successful realizations of instability models have concentrated on evolving the 
four giant planets onto their current orbits from a more compact configuration. Simultaneously, 
the possibility of forming and ejecting additional planets has been discussed, but never successfully 
implemented. Here we show that a large array of 5-planet (2 gas giants -I- 3 ice giants) multi-resonant 
initial states can lead to an adequate formation of the outer solar system, featuring an ejection of an 
ice giant during a phase of instability. Particularly, our simulations demonstrate that the eigenmodes 
which characterize the outer solar system's secular dynamics can be closely matched with a 5-planet 
model. Furthermore, provided that the ejection timescale of the extra planet is short, orbital excitation 
of a primordial cold classical Kuiper belt can also be avoided in this scenario. Thus the solar system is 
one of many possible outcomes of dynamical relaxation and can originate from a wide variety of initial 
states. This deems the construction of a unique model of solar system's early dynamical evolution 
impossible. 

Subject headings: Planets and satellites: dynamical evolution and stability; Kuiper belt: general 



1. INTRODUCTION 

In the last two dec ades, discoveries of the Kuiper belt 
(Jewitt & Luu 1993 ), as well a s planets orbiting stars 
other than the Sun (Mayor & Qucloz 1995), have sup- 
plied the centuries-old quest to understand the formation 
of the solar system with fresh constraints and insights 
into physical processes at play. Among a multitude of 
newly proposed format ion scenarios, the "Nice" model 
(iTsiganis et al. 2005; G omes et al.|2005[|Morbidem et al. 
2U05p is particularly notable, as it has attained a con- 



siderable amount of success in reproducing the various 
observed features of the solar system. Within the con- 
text of the scenario envisioned by the Nice model, giant 
planets start their post-nebular evolution in a compact, 
multi-resonant configuration, and following a brief period 



(Morbidelh et al.|2007 


Batygin & Brown|2010 


Walsh et 


al. 2011). 



to quantitatively reproduce the observed orbits of the gi 
ant planets, as well as their dynam ical architecture (i.e. 
secular eigenmodes of the system) (Tsiganis et al. 2005 
Morbidelli et al. 2009). Simultaneously, the brief msta- 
bility, inherent to the model, pro vides a natural trig ger 
to the Late Heavy Bombardment ( Gomes et al.|2005 ), as 
well as a transport mechanism for emplacement ot dy- 
namically " hot" Kuiper belt ob jects (KBOs) from inside 
of ^ 35 AU ( Levison et al.|2008 ). Meanwhile, it has been 
recently demonstrated that survival of a dynamically 
"cold" primordial population between Neptune's current 
3:2 and 2:1 exterior mean-motion resonances (MMRs) is 
fully consistent with a Nice model-like evolution of the 
planets, implying an in-situ formation of the cold classi- 



cal population of the Kuiper belt ( Batygin et al.j|2011 ). 
Finally, the presence of Jupiter's and Neptune's I'rojan 
asteroids has been attributed to chaotic capture of plan 



etesimals during the instability ( Morbidelli et al. 2005 



2007) 



Nesvorny et al. 

All successful realizations of the Nice model to date 
have been comprised exclusively of the four currently 
present giant planets. However, there exists no strong 
evidence that suggests that additional planets were not 
present in the solar system, at the epoch of the dispersion 
of the nebula. In fact, the oretical arguments, presented 
by Goldreich et al. (2004) point to a possibility of ini- 
tially forming as many as five ice giants, three of which 
get subsequently removed via ejections (see however Lev- 



ison fc Morbidelli (2007)). The dynamical sensibility of 
such a scenario is ifurther strengthened by the fact that a 
considerable fraction of standard Nice model simulations 
result in an ejection of an ice-giant after an encounter 
with at least one of the gas giants. 

In this paper, we explore an instability-driven dynam- 
ical evolution of a 5-planet system (2 gas giants -I- 3 ice 
giants) with an eye towards identifying a pathway to- 
wards reproduction of solar system-like dynamical archi- 
tecture. In principle, the realm of possibility available to 
this study is enormous. Consequently, rather than per- 
forming a comprehensive parameter-search, here we limit 
ourselves to systems that contain an additional Uranus- 
like planet, with the aim of presenting a few proof-of- 
concept numerical experiments. The plan of the paper is 
as follows. In section 2, we describe our numerical setup. 
In section 3, we show that in our model, the planetary 
orbits, their secular eigenmodes, as well as various popu- 
lations of the Kuiper belt are approximately reproduced. 
We conclude and discuss our results in section 4. 
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Fig. 1. — Orbital evolution of planets. Eaeh planet's semi-major axis, as well as perihelion and apohelion distances are shown as functions 
of time. The actual perihelion and apohelion distances of the planets are also shown for comparison as black error bars. In both cases, 
the innermost ice-giant ejects during the transient phase of instability, leaving behind 4 planets, whose orbits resemble that of the solar 
system. See main text for a description of the initial conditions. 



2. NUMERICAL EXPERIMENTS 

The numerical setup of the simulations perfor med here 
was quali t ativel y sin iilar to those present ed by IBatygm 



fc Brown] (|2010[) and [Batygin et all (|2011[) . Particularly, 



the hve giant planets were initialized m a compact, multi 
resonant initial condition, surrounded by a massive plan- 
etesimal disk that extended between its immediate sta- 
bility boundary and 30AU. 

Two of the three ice-giants were taken to have the same 
mass as Uranus and were initially placed next to Saturn 
and as the outermost planet respectively. The middle 
ice-giant was taken to have Neptune's mass. In all sim- 
ulations, Jupiter and Saturn started out in a 3:2 MMR, 
in accord with the results of hydrodynamical simulations 



ula (iMasset & Snellgrove|2001 


Morbidelli & Crida|2007 


^iercns & JNcison 20085. i'he ice giants were also se- 



ing dissipat ive forces, designed to mimic the presence of 
the nebula ( Lee fc Peale||2002 ) . Following resonant lock- 



ing, each assembled muiti-resonant initial condition was 
evolved in isolation for lOMyr, as an immediate test of 
orbital stability. This procedure yielded a total of 81 
stable multi-resonant initial conditions. 

The search for adequate dynamical evolutions was per- 
formed in two steps. First, we evolved 10 permutations of 
each initial condition, with planetesimal disks composed 
of = 1000 planetesimals. Disk masses were chosen 
randomly between = 25Me and M'^^^ = lOOM®. 

The density profiles followed a power-law distribution, 
E oc r*^ where the power-law index, fc, was chosen ran- 
domly between fcmin — 1 and fcmax = 2. The plan- 
etesimals were initialized on near-circular orbits (e ^ 
sin? ~ 10~^). To reduce the already substantial com- 
putational cost, self-gravity of the planetesimal swarm 
was neglected. 

Subsequently, we eliminated all initial conditions that 
did not yield any final systems that were comprised of 4 
planets, reducing the number of viable initial conditions 
to 25. Then, an additional 30 permutations of these ini- 
tial conditions were integrated with disks composed of 
N = 3000 planetesimals (but otherwise identical to those 



described above). Each integration was perfo rmed using 
the m ercury6 integration software package (Chambers 
19991 and spanned SOMyifM The calculations were per- 
tormed on Caltech's PANuU super-computer. 

After their completion, simulations that were deemed 
successful were reintegrated with the use of tracer siin u- 
lations (see 
to address 



Levison et al. ( |2008[ ); [Batygin et aL| ( |2011[ )), 
;he dynamical evolution of a locally formed 
population of KBOs. In particular, each run was sup- 
plemented with an additional disk of mass-less particles 
that resided in the cold classical region of the Kuiper 
belt (i.e. between the final exterior 3:2 and 2:1 MMRs of 
Neptune). 

3. RESULTS 

Out of the 810 integrations that were initially per- 
formed, 214 (~ 25%) cases featured an ejection of a single 
ice-giant, yielding a system composed of 4 planets. Per- 
haps unsurprisingly, in most cases the ejected planet is 
the ice-giant that neighbors Saturn. Of the 750 simu- 
lations that were performed following the elimination of 
initial conditions, 33 evolutions resulted in orbits remi- 
niscent of the solar system. Specifically, we searched for 
solutions where the Saturn- Jupiter period ratio exceed 
2, while the final semi-major axes of the ice-giants were 
within 3AU of their observed counterparts. No strong 
requirements were placed on the planetary eccentricities 
and inclinations. 

It is noteworthy that evolving the giant planets onto so- 
lar system-like orbits is insufficient for a simulation to be 
deemed successful for indeed, there are additional con- 
straints that must be satisfied. The first orbital con- 
straint is the reproduction of the secular architecture of 
the system. The secular orbital angular momentum ex- 
change (i.e. eccentricity evolution) of a planetary system 
containing N secondaries can be approximately repre- 
sented as a superposition of A^ eigenmodes, each cor- 
resp onding to a fundamental frequency of the syste m 
(see [Murray fc Dermott] ( [1999[ ); [Morbideni et alT] ( |2009[ )). 
Physically, the maximum eccentricity that a given planet 

^ Note that here, we make not attempt to time the onset of 
instability with the late heavy bombardment. 
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important constraint of the Nice model. [Levison et| 
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Fig. 2. — Eccentricity distribution of the remnant planetesimal disk. Red dots represent objects that have been dynamically emplaced, 
while the blue dots depict the locally formed cold classical belt at i = 50Myr. The red and blue triangles show objects whose orbits are 
stable on a 500Myr timescale. The scattered disk is shown with two solid curves and Neptune's exterior 3:2 and 2:1 MMR's are labeled with 
dashed lines. In the simulation presented in panel A of Figure (1), the cold belt suffers numerous encounters with the ejecting ice-giant, 
yielding considerable orbital excitation. Moreover, the inner cold belt is further dynamically depleted over SOOMyr of evolution. In the 
simulation presented in panel B of Figure (1), there is only a single close encounter between the cold belt and the ejecting ice-giant, yielding 
a dynamically cold orbital structure. 

attains in its secular cycle is equal to the sum of all of 
its corresponding eigenmode amplitudes. In the con- 
text of the traditional Nice model, some difficulty has 
been noted in correctly reproducing the dynamical char- 
acter of Jupiter's and Saturn's eccentricity evolution. 
Particularly, it has been shown that smooth passage of 
Jupiter and Saturn through the 2:1 MMR has a ten- 
dency to under-excite the ec centricity eigenmode a s 
well as their mutual inclinations ( Morbidelli at al.|2009 ). 
This difficulty can be overcome sometimes by invoking a 
close-encounter between an ice-giant and the gas giants. 
Consequently, we have checked, u s ing F ourier analysis 



(see for example Morbidelli et al. (|2009[)), the relative 



strength of the and eigenmodes m all simulations 
whose orbital end-state resembled the solar system. We 
did not restrict the success criteria of our simulations to 
include the correct reproduction of the mean eccentric- 
ities of the planets and in some cases, the mean final 
eccentricities (namely those of the ice-giants) exceeded 
their observed counter-parts by as much as a factor of 
~ 2. This is, however, likely an artifact of the coarse 
representation of the planetesimal disk and the resulting 
dynamical friction, that we employed in our calculations 
and should not be viewed as a major drawback. 

In total, we found ten cases (corresponding to eight dif- 
ferent initial conditions) where their amplitudes are sat- 
isfactorily reproduced. Specifically, in these ten cases, 
the amplitude of Jupiter's g^ eigenmode exceeds that 
of the (76 eigenmode, while the amplitudes are roughly 
the same for Saturn. We made no attempt at quantita- 
tively matching the pair's inclination eigenmodes, how- 
ever their repro duction does not appear to be problem- 
atic ( Morbidelli et al.||200 9). We also examined the am- 
plitudes of the secular eigenmodes of Uranus and Nep- 
tune. Generally it appears that the dynamical archi- 
tecture of the ice-giants is set in an essentially random 
manner, depending on the particular encounter history. 
Consequently, we decided to not use ice-giant secular ar- 
chitecture as a distinctive property in our analysis. 

Successful formation of the Kuiper belt is another 



al. (2008) have shown that the excited populations of 
the Kuiper belt are naturally emplaced from inside of 
^ 35AU during the instability. Given that this aspect 
of the planetary evolution is not particularly different 
between the 4-planet and the 5-planet scenarios, there 
is little reason to speculate that the dynamical pathway 
for formation of the resonant, scattered, and hot classi- 
cal populations of the Kuiper belt will be inhibited. The 
same is likely to be true for chaotic capture of Trojan 
asteroids. The cold classical population of the Kuiper 
Belt, however, is a different story. 

A series of observational dissimilarities between the 
cold classical population and the rest of the Kuiper 



belt (e.g. uniq uely red colors (Trujillo & Brown 



2002 



Lykawka fc Mu kai 200 51); strongly enhanced wide bi 
nary fraction ( Stephens & Noll|2006 Parker & Kavelaars 
2010)) suggest that the cold classicais formed in situ, and 
maintained dynamical coherence despite Neptune's tem- 
porary acquisition of high eccentricity and inclination (a 
characteristic orbital feature of the cold classica l popu- 
lation is inclination that does not exceed ~ 5° (Brown 



2001)). In a recent study, Batygin et al. (2011) showed 
that local formation of the cold classicais is lully consis- 
tent with an instability-driven evolution of the planets, 
given favorabl e conditions d uring the instability. Par- 
ticularly, Batygin et al. ( 2011J required the apsidal pre- 
cession and nodal recession rates of Neptune to be com- 
paratively fast to prevent secular excitation of the cold 
classical orbits, in addition to a sufficiently small apohe- 
lion distance of Neptune, to avoid orbital excitation due 
to close encounters. 

In a 5-planet scenario, the retention of unexcited orbits 
of the cold population can be jeopardized by the ejecting 
planet. This possibility served as a premise for recalcu- 
lation of the dynamical evolution of test particles in the 
cold classical region, with the aid of tracer simulations. 
We found that in only three of the ten simulations where 
the secular eigenmodes were successfully reproduced, a 
primordial cold classical population of the Kuiper belt 
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Fig. 3. — Inclination distribution of the remnant planctesimal disk. Red dots represent objects that have been dynamically emplaced, 
while the blue dots depict the locally formed cold classical belt at i = 50Myr. The red and blue triangles show objects whose orbits are 
stable on a 500Myr timescale. The tentative i = 5° boundary between the cold and hot classical belts is shown with a solid line, while 
Neptune's exterior 3:2 and 2:1 MMR's are labeled with dashed lines. In the simulation presented in panel A of Figure (1), the cold belt 
suffers numerous encounters with the ejecting ice-giant, yielding considerable orbital excitation. Moreover, the inner cold belt is further 
dynamically depleted over SOOMyr of evolution. In the simulation presented in panel B of Figure (1), there is only a single close encounter 
between the cold belt and the ejecting ice-giant, yielding a dynamically cold orbital structure. 



also retained unexcited orbits. Consequently, it appears 
that the retention of an unexcited cold belt is not the 
norm of the 5-planet scenario. This is not surprising, 
since dynamical excitation can only be avoided if ejec- 
tion and close encounter timescales are sufficiently short. 
From a quantitative point of view, our simulations imply 
that dynamically cold orbits can only be sustained if the 
ejecting ice-giant spends ~ 10'* years or less crossing the 
classical Kuiper belt region. As a result, we can expect 
that the incorporation of yet additional ice- giants into 
the p rimordial solar system (as suggested by Goldreich 
et al. ( j2004^ ) will further diminish the chances ot repro- 
ducing the cold classical Kuiper belt, since all ejecting 
planets would have to do so very rapidly. 

Two of the three successful simulations discussed above 
are presented in Figures (1-3). Specifically, Figure (1) 
shows the orbital evolutions of the runs, while Figures 

(2) and (3) show the eccentricity and inclination distri- 
butions of planetesimals in the Kuiper belt region respec- 
tively. The red and blue dots, shown in Figures (2) and 

(3) depict the orbital distribution of emplaced and local 
planetesimal populations respectively at t = 50Myr. The 
red and blue triangles depict the planetesimals whose or- 
bits remain stable at i = 500Myr. 

The starting multi-resonant initial condition of the 
simulation presented in panels (A) of the figures is one 
where Saturn and the first ice-giant are locked in a 2:1 
MMR, while both pairs of ice-giants are locked in 5:4 
MMRs. The planetesimal disk in this simulation was 
comprised of iV = 3000 particles and contained a total 
of 26M0. In this evolution, the local population of test 
particles suffers numerous short close encounters with the 
ejecting ice giant, yielding a more excited and depleted 
cold classical population, compared to that of run (B). 
Note also that at t = 50Myr, the ice-giant eccentrici- 
ties are considerably greater than that of Uranus and 
Neptune. These high eccentricities do not get damped 
away by dynamical friction in the following 500Myr of 
dormant evolution. Consequently, in this simulation, the 
inner edge of the cold belt gets dynamically depleted over 



the following 500Myr. It is furthermore noteworthy that 
another simulation that originated from the same initial 
condition reproduced the eigenmodes of the system cor- 
rectly, although the primordial cold Kuiper belt in this 
integration was entirely destroyed by close encounters. 

The starting multi-resonant initial condition of the 
simulation presented in panels (B) of the figures is one 
where Saturn and the first ice-giant are locked in a 3:2 
MMR, the inner pair of ice-giants is locked in a 2:1 
MMR, and the outer pair of the ice-giants is locked 
in a 4:3 MMR. In this simulation, the disk consisted 
of iV = 1000 particles and had a cumulative mass of 
42M0. Incidentally, the frequency spectrum of the ec- 
centricity vectors of Jupiter and Saturn produced in this 
simulation, matches that of the real Jupiter and Saturn 
exceptionally well, signaling a nearly ideal reproduction 
of the secular eigenmodes. Particularly, the simulatio n 
yields (in the notation of (Murray & Dermott 1999)) 

(sim) / (sim) 



'55 



/e. 



56 



2.28 and 655 /ck 



'66 



55 



0.0465, el^g" = 0.06 7 while the 
0.0442, 666 = 0.04 82 (|Morbidelli 



0.51 where as 

the solar system is characterized by 655/656 = 2.81 and 
665/666 — 0.68. The scaling of the eigen- vectors are also 
well reproduced: e 
solar system has 655 

et al. 20091. In this simulation, the cold liuiper belt 
sutters only a single short encounter with the escaping 
ice-giant, allowing for the orbits (inclinations in particu- 
lar) to remain dynamically cold. 

4. DISCUSSION 

In this paper, we have presented a successful realiza- 
tion of the Nice model which starts out with 5 plan- 
ets. The numerical experiments presented here explic- 
itly show that such an evolution is plausible since the 
resulting 4-planet systems can closely resemble the so- 
lar system. Particularly, in both simulations presented 
here, the secular architecture of the outer solar system is 
well reproduced. Furthermore the demonstrated survival 
of a local, primordial cold classical Kuiper belt suggests 
that all constraints that can be matched with a 4-planet 
model can also be matched with a 5-planet model to an 
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equal degree of satisfaction. 

It is noteworthy that ejection was not always necessary 
in our simulations to generate a 4 planet system. In a 
handful of runs (one of which successfully reproduced the 
secular eigenmodes, but not the cold Kuipcr belt), one of 
the ice-giants ended up merging with Saturn. In princi- 
ple, such a scenario may help explain Saturn's enhanced 
metallicity in comparison with Jupiter . Although, here 
again the explanation is not unique (see Stevenson ( 1982 1 
and the references therein). 

In a traditional realization of the Nice model, the rate 
of successful reproduction of the secular eigenmodes is 
rather low i.e. ~ 10 % of the integrations for a favor- 
able initial condition (Batygin & Brown 20101. This is 
in part because an ice-giant/gas giant encounter often 
leads to an ejection of the ice-giant, leaving behind only 
three planets. Thus, the need for an ice-giant/gas-giant 
encounter in the orbital history of the solar system is in 
itself motivation for a 5-planet model. 

The statistics of simulations presented in this work 
suggest that a 5-planet model is neither more nor less 
advantageous. Recall that the probability of ending up 
with only 4 planets is 214/810 25%. The probabil- 
ity of reproducing the secular eigenmodes of Jupiter and 
Saturn is 10/750 ~ 1.5%. Naively, this yields an overall 
probability of success of only ~ 0.4%. However, it is im- 
portant to keep in mind that the characteristic outcomes 
are generally dependent on initial conditions]^ and runs 
that originated from the initial condition presented in 
panels (A) correctly reproduced the secular eigenmodes 
in 2 out of 30 simulations. This statistic is similar to the 

4- planet model. Finally, the 1/30 probability of also re- 
taining an unexcited cold classical Kuiper belt puts the 

5- planet model and the 4-planet model on equal footing 



in terms of success rate ( Batygin et al.|20"TT |. That said, 
it is important to note that this success rate is only char- 
acteristic of the particular 5-planet model that we have 
constructed. In other words, it is likely that if one allows 
the mass of the ejected planet to also be a variable pa- 
rameter, tuning of the initial state may in principle lead 
to a more frequent reproduction of the solar system. 

The results presented in this work imply that the solar 
system is one of many possible outcomes of dynamical 
evolution, and can originate from many possible initial 
conditions. As a result, the possibility of having an extra 
planet initially present in the system, yet its ejection 
leaving no observable signature erases any hope for 
construction of a deterministic model for solar system 
evolution. The forward process-like nature of the Nice 
model is not surprising, given that the solar system 
exhibits large-scale chaos, characterized by Lypunov 
times that are comparable to orbital timescales, during 
the instability. Moreover, the similarity between the 
orbital architectures of simulations whose outcomes 
were deemed unsuccessful in this work and those of 
extra-solar planetary systems further confirms that 
planet-planet scattering is likely to be the physical 
process res ponsible for shaping the orbital distribution 
of planets ( Juric & Tremaine 2008 ) . Consequently, we 
conclude that an instability-driven dynamical history 
remains a sensible choice as a baseline scenario for solar 
system's early dynamical evolution. 
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